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Abstract. The eruption of V4332 Sgr discovered in February 1994 shows striking similarities to that of V838 Mon started in 
January 2002. The nature of these eruptions is, however, enigmatic and unclear. We present new photometric and spectroscopic 
data on V4332 Sgr obtained in April-May 2003 at the SAAO. The obtained spectrum shows an unusual emission-line compo- 
nent superimposed on an early M-type stellar spectrum. The emission-line spectrum is of very low excitation and is dominated 
by lines from neutral elemets (Nal, Fel, Cal) and molecular bands (TiO, ScO, AlO). We also analyse all the observational data, 
mainly photometric measurements, available for V4332 Sgr This allows us to follow the evolution of the effective temperature, 
radius and luminosity of the object since February 1994 till 2003. We show that the observed decline of V4332 Sgr can be 
accounted for by a gravitational contraction of an inflated stellar envelope. The combined optical and infrared photometry in 
2003 shows that apart from the M-type stellar component there is a strong infrared excess in the KLM bands. This excess 
was absent in the 2MASS measurements done in 1998 but was probably starting to appear in K in 1999 when the object was 
observed in the DENIS survey. We interpret the results in terms of a stellar merger scenario proposed by Soker & Tylenda. The 
infrared excess is likely to be due to a disc-like structure which is either of protostellar nature or has been produced during the 
1994 eruption and stores angular momentum from the merger event. 
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1. Introduction 

V4332 Sagittarii was discovered towards the end of 
February 1994 as an apparent nova. Subsequent spectroscopic 
observations however indicated that this was not a classical 
nova. It lacked spectral features characteristic of novae while 
its spectral type, K3^ III-I on March 4, was too late for a 
nova at maximum and with time, contrary to classical novae, it 
evolved toward lower and lower effective temperatures. A de- 
tailed description and analysis of available observational data 
for V4332 Sgr in eruption can be found in Martini et al. (1999, 
hereafter MWT99 1. 

Similar evolution has been observed for V838 Mon which 
was discovered in eruption in early January 2002. This object 
underwent more complex photometric and spectral evolution, 
but in March- April 2002 its behaviour was quite similar to that 
observed for V4332 Sgr (e.g. Munari et al. 2002 Crause et al. 
l2003l . Therefore V838 Mon and V4332 Sgr (and presumably 
also M31 RV - see Mould et al.'1990) are considered a myste- 
rious new class of eruptive variables. As discussed in Soker & 
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Tylenda ('2003'), thermonuclear events, i.e. nova-type runaway 
or late He-shell flash, cannot account for the observed evolu- 
tion of these objects. Instead these authors propose that the 
eruptions can be due to stellar merger events (see also Retter 
& Marom 2003 1. 

V4332 Sgr has been forgotten by observers since its decline 
in June 1994. As far as we know, the only observations obtained 
between 1994-2003 are 2MASS measurements in May 1998 
and DENIS photometry in September 1999. Motivated by the 
interest in V838 Mon, we took low resolution optical spectra 
of V4332 Sgr in April 2003, followed by photomertic measure- 
ments in May 2003. We were surprised to find that the object 
displayed an unusual emission-line spectrum, completely dif- 
ferent to that of typical emission-line objects such as planetary 
nebulae, symbiotic stars or T Tauri stars. The strongest emis- 
sion feature appeared to belong to the Na I D Unes, but later we 
realised that many emission lines are due to molecular bands of 
TiO, AlO and ScO. Banerjee et al. (2003, hereafte r BVA03 ) ob- 
served V4332 Sgr in near-IR in April-June 2003. Their spectra 
also show strong emission features which they attribute to AlO. 
More recently, in September 2003, Banerjee & Ashok (2004, 
hereafter lBA04t observed the object in the optical and obtained 
a spectrum very similar to ours. Very recently, when the orig- 
inal version of this paper was already submitted for publica- 
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Table 1. Cousins photometry of V4332 Sgr obtained on 21 
May 2003 



UT 


Band 


Magnitude 


Error 


03:50 


B 


19.6 


±0.2 


03:15 


V 


17.63 


±0.03 


03:27 


Rc 


16.34 


±0.00 


03:39 


h 


15.09 


±0.00 



tion, Banerjee et al (2004, hereafter lBVA04t presented infrared 
spectra obtained in September 2003 and April 2004. They have 
detected a water-ice absorption band and a CO band in emis- 
sion. Thus V4332 Sgr is now an unusual emission-line object 
that deserves detailed observations. In this paper we present our 
observations and analyse the available data on the evolution of 
V4332 Sgr since its eruption in 1994. 



2. Observations 

The observations of V4332 Sgr were made at the South African 
Astronomical Observatory (SAAO). The spectra were obtained 
on 3 April 2003 with the grating spectrograph with SlTe CCD 
mounted at the Cassegrain focus, f/18, of the 1.9-meter tele- 
scope. The total duration of the exposure was 900 seconds with 
the spectral coverage 3500-7400 A and an average resolution 
of 1000 using grating No. 7 (300 lines/mm). The selected slit 
width was 300jum giving I'.'S projected on sky, comparable to 
seeing conditions during the observations. The FWHM of the 
night sky lines was 5.6 A. Two standard stars, CD-32°9927 
and EG 274, were observed during the night. This allowed us 
to calibrate the spectrum in flux units while the wavelength cal- 
ibration was secured by observing CuAr lamp spectra. 

The standard reduction procedures were applied to the 
spectra with the MIDAS' longslit package. These included bias 
subtraction, flat-field correction, atmospheric extinction cor- 
rection, wavelength and flux calibration. 

Cousins BVRcIc photometry was obtained with the 1.0- 
meter telescope and the Ik x Ik STE4 CCD camera on 
21 May 2003. The frames were bias-subtracted, overscan- 
trimmed and flat-fielded using the CCDPROC package within 
IRAF. Target observations were calibrated against Landoldt 
standard stars. The results are given in Tabled Note that the er- 
rors given in the table are based only on the photon noise, they 
do not take into account external sources such as the calibra- 
tion. Thus the true errors are probably larger than those quoted 
in the table. 

The results in Table [l] can be compared to those obtained 
bv lBA04l on 29 September 2003. Apart from the B magnitude, 
which is uncertain in both measurements and differs by 0.4, the 
other magnitudes are consistent within 0.03 - 0.11. Thus we 
can conclude that the object did not evolve significantly during 
the 4 month period in 2003. 



' MIDAS is developed and maintained by the European Southern 
Observatory 



3. Spectrum 

The resultant spectrum of V4332 Sgr obtained on 3 April 2003 
is presented in Fig.^ The observational noise becomes exces- 
sive shortward of ~4500 A where the spectrum is faint so that 
section is not shown in Fig.E For the remaining part a contin- 
uum level has been fitted by the median filtering of the spec- 
trum and the observational noise has been derived by calculat- 
ing the rms of the deviations from this continuum, excluding 
the apparent emission features. The calculated average noise 
is 0.19 X 10"'^ erg cm"^ s"' A"' and as a result the signal to 
noise ratio (S/N) for the detected continuum changes from only 
1 to ~5 at the longest wavelenghs. However for the emission 
lines the S/N can be substantially higher. Our spectrum can be 
compared to that of BA04 Note, however, that our wavelength 
range does not extend far enough into the red to register the 
prominent lines of KI. 

As Fig. n shows, the spectrum is dominated by prominent 
emission features. This emission-line spectrum is, however, un- 
usual and does not resemble any typical emission-line objects, 
such as planetary nebulae, symbiotic stars or T Tauri stars. Its 
is also significantly different from that displayed by V4332 Sgr 
during the 1994 eruption <MWT99t when the emission-line 
spectrum was dominated by the Balmer series. To some extent 
the present spectrum of V4332 Sgr can be compared to those 
of pecular giants, such as U Equ (e.g. Barnbaum et al. ll996t or 
VY CMa (e.g. Herbig [T974l . 

The strongest emission features in the spectrum of 
V4332 Sgr are identified in Fig. [0 Table |2l presents details 
of the emission lines, i.e. measured wavelengths (in A), line 
fluxes (in 10"'^ erg cm"^ s"') with errors in parenthesis, 
full width half maxima (in A), equivalent widths (in A) and 
the proposed identifications. For all lines the flux, FWHM 
and EW were derived by performing a Gaussian fit to the 
observed profile using the REWIA package (J. Borkowski; 
http://www.ncac.torun.pl/~juborki. The errors on the line fluxes 
mostly depend on the uncertain continuum level and have 
therefore been calculated in an aproximate way as a product 
of the half of the average noise level and a doubled FWHM of 
a given line. 

The strongest emission feature in our spectrum of 
V4332 Sgr is due to the unresolved Nal D lines. We have also 
been able to identify emission lines from other neutral atoms, 
i.e. Fe I and Cal (see below). 

The most interesting aspect of the present spectrum of 
V4332 Sgr is that it displays molecular bands in emission. The 
strongest and most numerous features are due to TiO. We have 
identified the prominent emissions near 7 1 00 A with three com- 
ponents of band (0-0) of the -y system corresponding to the tran- 
sitions between electronic levels A - X^A. We believe that 
the (1-1) band of the same system at 7 100-7200 A (not marked 
in Fig. ^ is also visible in our spectrum. Another prominent 
band of TiO is that of (0-0) of the y' system (B^U - X ^A) 
near 6200 A. Some other emission features can also be identi- 
fied with excited bands of the TiO y system. Following Herbig 
(I1974I I we have identified the prominent emission features at 
6040 and 6080 A with ScO. AlO, which dominates the near-IR 
spectrum of .BVA03. is also present in emission in our spec- 
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Fig. 1. The spectrum of V4332 Sgr obtained on 3 April 2003. The observed flux (in 10"'^ erg cm"^ s ' A"') is shown versus the 
wavelength (in A). Due to the excessive observational noise, the spectrum below ~4500 A is not shown. The sharp spike near 
5580 A is due to imperfect sky subtraction. The labels identify the strongest emission features. 



tram. Four bands of the B - X electronic transitions are 
cleai-ly seen near 4850 A and 5100 A. 

The strong feature at 6569.8 A was initially attributed to 
H„, but we subsequently identified it as the Ca 1 intercombina- 
tion transition. H„ was seen in the spectrum of V4332 Sgr dur- 
ing its 1994 eruption (MWT99l, but at that time its mean posi- 
tion was 6558.7 A. Since our position of the Nal D lines differs 
only by 1.3 A from that of MWT99I it is difficult to reconcile 
a difference of 1 1 A in the Ha position. |^^04 suggest that this 
feature might be due to TiO y' (0,1) as the wavelength matches 
well. In that case, other emissions originating from the same 
transition should be seen at 6596.3 A and 6629.0 A but they are 
not present. Furthermore, Table|2]shows that the 6569.8 A fea- 
ture is significantly narrower than other molecular lines, sug- 
gesting that it is an atomic line. Given the very low excita- 
tion of the spectrum of V4332 Sgr, we consider Cal 6572.8 A 
('S-^P") as the most probable identification for the 6569.8 A 
feature. 

As can be seen from Fig. ^ the spectrum shows a clear 
continuum rising toward longer wavelengths. The rather low 
signal-to-noise ratio does not allow unambiguous identification 
of absorption features. Yet a comparison of the observed spec- 
trum with synthetic ones suggest an early M type. 

An analysis of the strengths and widths of molecular bands 
allows us to estimate excitation conditions in the medium 
where the bands are produced. For the TiO bands we have used 
the data base described by Schwenke ( 1998 ) and distributed by 
Kuracz <1999t . The observed widths of the emission peaks of 
the (0-0) band of the y system of TiO (7050-7 130 A) has al- 
lowed us to estimate the rotational temperature, Ti-ot - 200 K. 
On the other hand, the relative equivalent widths of these emis- 
sion peaks point to the electronic temperature, Tei ^ 500 K. 
Following the analysis of Herbig Cl974j and using the ob- 



served equivalent- width ratio of the 6079 A and 6036 A peaks 
we can estimate T^ot - 800 K for ScO. It is worth noting that 
from numerous bands of AlO in the near-IR l5VA03l estimated 
Tvib - 3000 K and T,ot - 300 K. On die other hand, from the 
water ice feature at 3.05 fj.m BVA04l have estimated a temper- 
ature of 30-50 K, while their CO band analysis has suggested 
300-400 K. 



4. Analysis of the photometric data 

Using the available photometric data we have attempted to 
analyse the evolution of V4332 Sgr since its outburst in 1994. 



The method is the same as that used in Tylenda i2005t in the 
analysis of V838 Mon. Basically it fits standard colours from 
a sequence of spectral types to the photometric data. From the 
best fit, obtained using the least squares method, one derives 
the effective temperature, Teff , and the angular stellar radius, 6. 
This procedure requires the value of interstellar extinction to 
be known. Finally, if the distance to the object is known, the 
linear stellar radius, R, and the stellar luminosity, L, can then 
be calculated. 

4.1. Interstellar extinction 

The interstellar extinction to V4332 Sgr has been determined 
bv lMWT99 Comparing the photometric colours with the spec- 
tral classification from their spectra obtained in March 1994 
they have found Eb-v = 0.32 + 0.02. However, the values 
of Ey-R and Ey-i they obtained were significantly bluer, even 
negative in most instances (see thek Fig. 6). This troublesome, 
at first sight, inconsistency can, however, be easily explained. 
The intrinsic colours which lM\VT99l compared with the photo- 
metric observations were in the Johnson U966j system while 
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Table 2. Emission line spectrum of V4332 Sgr observed on 
3 April 2003. [^obs, FWHM and EW are in A. Flux is in 
10^'^ erg cm"^ s"'. Identification gives A].^\, (A), element 
(molecule) and transition (for molecular transitions (v'-v") 
means the upper and lower vibrational levels).] 



'lobs 


Flux 


FWHM 


EW 


Identification 


4606.7 


7.6 (+3.3) 


17 


82 




4846.9 


7.8 (±2.0) 


11 


52 


4842.3 AlO -X-S* (0-0) 


4866.7 


6.8 (±2.4) 


12 


44 


4866.4 AlO -X'!,^ (l-l) 


5078.1 


8.4 (±2.3) 


12 


41 


5079.4 AlO S -S+ - ;C -E+ (0- 1 ) 


5105.9 


9.4 (±2.6) 


14 


44 


5102.1 AlO (1-2) 










5110.4 Fel (1) 


5123.1 


3.3 (±1.7) 


9.2 


15 


5123.3 AlO -X^I,^ (2-3) 


5164.3 


4.8 (±1.4) 


7.4 


21 


5166.3 Fel (1) 










5168.9 Fel (1) 


5204.3 


10.8 (±1.7) 


9.2 


46 


5204.6 Fel (1) 


5294.4 


3.8 (±1.2) 


6.3 


15 




5802.2 


7.4 (±3.9) 


21 


19 




5891.0 


64.5 (±2.2) 


12 


155 


5889.9 Nal 










5895.9 Nal 


6039.6 


12.9 (±3.0) 


16 


29 


6036. 1 ScO A-n- X L (0-0) 


6079.5 


14.2 (±3.3) 


17 


31 


6079.2 ScO A^TI-X^I. (0-0) 


6137.3 


5.5 (±3.8) 


20 


12 




6161.2 


12.4 (±2.9) 


15 


26 


6161.6 TiO y (0-0) 


6190.7 


6.5 (±2.4) 


13 


14 


6185.6 TiO y (0-0) 


6221.6 


8.3 (±3.6) 


19 


17 


6216.8 TiO y (0-0) 


6403.5 


9.0 (±4.1) 


22 


18 




6569.8 


15.7 (±1.6) 


8.4 


32 


6572.8 CaU'S-'P") 


6780.3 


6.8 (±3.4) 


18 


14 


6781.3 TiO 7 (2-1) (?) 


6840.3 


9.3 (±6.5) 


34 


19 


6849.9 TiO 7 (3-2) (?) 










6814.7 TiO 7 (3-2) (?) 


6974.8 


10.2 (±4.0) 


21 


21 




7057.5 


26.3 (±3.8) 


20 


52 


7054.3 TiO 7 (0-0) 


7093.5 


24.1 (±3.3) 


17 


47 


7087.7 TiO j (0-0) 


7129.4 


30.4 (±2.7) 


14 


57 


7125.5 TiO 7 (0-0) 










7124.9 TiO 7 (1-1) 


7160.9 


4.7 (±2.2) 


11 


8.6 7158.5 TiO 7 (1-1) 


7199.4 


5.5 (±2.4) 


13 


9.5 7197.2 TiO 7 (1-1) 










7197.2 TiO 7 (2-2) 


7271.2 


3.2 (±3.0) 


16 


4.8 


7270.4 TiO 7 (2-2) 


7343.2 


5.7 (±2.9) 


16 


6.6 





the observations of Gilmore were done in the Cousins (R^ and 
Ic) system (Gilmore, private communication). Indeed, as dis- 
cussed in Sect. l4.3l and l4.4l there is a good agreement between 
the spectral types of MWT99 and the spectral types resulting 
from our fits to the BVRJc magnitudes if Eb-v - 0.32 is as- 
sumed. We can also conclude that the uncertainty of this value 
is not greater than 0. 1 . For Eg-v beyond this range our spectral 
classes from the photomertic data become systematically dis- 
crepant from those of lMWT99l In particular, from the fits done 
to the photometric data on 4-5 March (these data have been 
used by MWT99 to compare with their spectrum on 4.5 March 
and to derive Eg-v = 0.32) assuming Eb-v = 0.20 we have ob- 
tained K4.5-6.5 while for Eb-v = 0.45 the result is K2.0-3.0. 
This can be compared to K3-4 obtained by MWT99.fi"om their 
spectrum on 4.5 March. 

In the present paper we assume Eb-v - 0.32. The uncer- 
tainty in Eb-v introduces uncertainties in the efi"ective temper- 



atures and the angular stellar radius derived from the fits dis- 
cussed below. This errors can next propagate to the values of 
the linear stellar radius and luminosity. We have found that an 
error of 0.1 in Eb-v afifects the eff'ective temperature and the 
angular stellar radius by at most 5%. Thus in most cases this 
source of uncertainties is small compared to that due to the 
reference spectra (giants or supergiants) used in the fitting pro- 
cedure (see Table|3}. The resulting error in luminosity is larger, 
i.e. 15-20%, but even larger uncertainties are expected to result 
from the distance uncertainty (see Sect. 14. 2t . 

4.2. Progenitor and distance 

MWT99 have estimated that the distance to V4332 Sgr might 
be ~300 pc assuming that the object was a K giant at maximum. 
They have however concluded that this procedure did not lead 
them to consistent results if later observations are considered. 

Skiff" J1994> . using the POSS prints, identified a possible 
progenitor of V4332 Sgr with B = 18 and = 16. Using the 
ALADIN service we have searched for brightness estimates of 
the V4332 Sgr progenitor in archival surveys. We have found 
data in three catalogues. The SuperCOSMOS catalogue for a 
star at the position of V4332 Sgr gives B = 18.66, R = 16.72 
and / = 15.89. In the USNO-Bl.O catalogue the figures are 
B = 18.20, R = 16.48 and / = 15.69.^ Finally in the USNO- 
A2.0 catalogue one finds B = 17.9 and R - 16.4. The results 
are fairly consistent although there are systematic difi'erences 
between the values in the three catalogues. In all the bands the 
object is faintest in SuperCOSMOS while brightest in USNO- 
A2.0. The colours are more consistent. This is not surprising 
as in authomatic measurements from survey plates or frames 
systematic errors are usually most important. 

With Eb-v - 0.32 (see Sect. 14. H the colours from the 
above measurements imply an effective temperature of 5000- 
6000 K. For this temperature range, the most reasonable as- 
sumption is that V4332 Sgr was a main sequence star prior 
to the eruption. Adopting intrinsic colours and absolute mag- 
nitudes for main sequence from Drilling & Landoldt ( 2000.) 
as well as assuming Eb-v = 0.32 and the standard extinction 
curve the magnitudes from SuperCOSMOS can be best fitted 
with a KO V star at a distance of 1.35 kpc. The same procedure 
applied to the results from USNO-B 1 .0 gives a G6 V type and a 
distance of 1.58 kpc while using the magnitudes from USNO- 
A2.0 one gets F8 V and 2.6 kpc. 

If one assumes that before eruption V4332 Sgr was a giant 
rather than a main sequence star, then the observed colours of 
the progenitor would imply a spectral type of G5-F8 III and 
Mv - +1.0 (Schmidt-Kaler 1982). Then the distance would 
increase to ~10-13 kpc. This solution seems, however, to be 
unlikely given the position of V4332 Sgr in the sky, I = 13'?63 
b = -9°.4Q, and the observed reddening, Eb-v - 0.32. Several 

2 Both, SuperCOSMOS and USNO-B 1 .0 give second R magnitudes 
for the V4332 Sgr progenitor, which are 14.96 and 14.38, respectively. 
From the R magnitudes of three nearby stars in these catalogues which 
in the red images are of similar brightness as the V4332 Sgr progeni- 
tor it can be safely concluded that these second magnitudes are signif- 
icantly too bright. 
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planetary nebulae near this position, i.e. PNG 11. 3-09.4, 12.5- 
09.8, 13.8-07.9, 13.7-10.6, 14.2-07.3, have Eb-v in the range 
0.37-0.52 (Tylenda et al. 1992), yet their distances are esti- 
mated to be within 3-8.5 kpc (Acker et al. .l992i) . 

In the present paper we assume that V4332 Sgr was a solar 
type star before eruption and that the object is at a distance of 
~1.8 kpc (mean value from the above estimates adopting main 
sequence). It is however clear from the above that this conclu- 
sion is subject to significant uncertainties. The distance value 
is probably uncertain to 50% or so introducing corresponding 
uncertainties in the values of the stellar radius and luminosity 
derived in Sect. 14.41 

4.3. Observational data and fitting them with standard 
spectra 

During the 1994 eruption, the only multicolour {UEVRch) 
photometry of V4332 Sgr was done by Gilmore (I 1994t and 
the data cover the period between 27 February and 1 1 March. 
[MWT99 have estimated the BVR magnitudes from their spec- 
tra obtained on 11 March and 5-6 June 1994. From the pub- 
lished spectra in MWT99 we have made rough estimates of 
the / magnitude. Their Fig. 4 shows that on 11 March the 
flux in the / band was ~3.5 times stronger than in the R 
band which results in R - I ^ 2.2. From Fig. 5 in MWT99 
one can estimate that on 5-6 June the flux in the / band was 
Fi ^ 3.5 X 10"'** ergs cm^^ s"' A"' giving / ^ 10.9. 

In the 2MASS survey the object was observed on 
18 May 1998 and flie results were: J = 12.10, H = 11.60 and 
K - 10.99. V4332 Sgr was also measured in the DENIS survey 
on 11 September 1999 which gave Icmm = 14.37, J = 12.46 
and K, = 10.66. 

In 2003, apart from our BVRJc photometry obtained on 
21 May and summarized in Tabled optical photometry was 
also done by BA^on 29 September. Besides, JHK magni- 
tudes were obtained by BVA03 on 19 June 2003, while LM 
magnitudes measured on 5 September 2003 can be found in 
[BVA04 As discussed in Sect. |2l V4332 Sgr did not evolve 
significantly between May and September Therefore we have 
combined all the data from 2003 into one set for the analy- 
sis discussed below. In particular, in the optical we have taken 
mean values from our results in Tableland those of lBA04l 

When fitting the observed magnitudes, we applied the stan- 
dard intrinsic colours for the luminosity class III (giants) and I 
(supergiants). This is justified for the period of the 1994 erup- 
tion as the observed spectra of V4332 Sgr then resembled those 
of the III-I luminosity class ( MWT99 ). Also, for the data ob- 
tained in 2003 we have found that the giant or supergiant spec- 
tra better fit the observations than the colours for main sequence 
stars. The intrinsic photometric colours have been taken from 
Schmidt-Kaler (fT9^ - UBV), Johnson ( fT966l - VRI for G and 
K types), Lee fWJU- VRI for M types), Bessell & Brett CM^ 
- VIJHKLM for giants) and Koornneef ( 1983,- VJHKLM 
for supergiants). The calibration of the effective temperature 
and the bolometric correction against the spectral type has been 
taken from Schmidt-Kaler (il982i. 




-0.4 -0.2 

log X(/j.m) 



Fig. 2. Evolution of the spectrum (magnitude versus effective 
wavelength) ofV4332Sgrin 1994. Symbols show the observed 
magnitudes (taken from Gilmore 1994 and MWT99 1. Solid and 
dotted curves show the fitted standard supergiant and giant 
spectra, respectively (reddened using the standard extinction 
curve and Eb-v - 0.32). From top to bottom the curves cor- 
respond to the observations made on 27.7 February, 4.7 March, 
7.7 March, 10.7 Mai'ch, 1 1 .4 Mai'ch and 5-6 June 1994, respec- 
tively. 



Figure 121 presents the best fits of the standard giant (dotted 
curves) and supergiant (solid curves) spectra to the photometric 
data obtained in 1994. Note that for clarity of the figure the 
spectra on 3.7 March and 8.7 March have not been displayed 
in Fig. 121 

The results of the fits, i.e. spectral types, effective tem- 
peratures and angular stellar radii, 6 (in radians), are given in 
Table |3l For each date the results of the fits are given in two 
lines: first the results from fitting the supergiant (I) spectra, sec- 
ond those for - the giants (III). The lines marked with M in the 
fourth column give the results from lMWT99 l i.e. spectral type 
from classification and effective temperature from stellar atmo- 
sphere model analysis of their spectra. 

Note that the U magnitude, measured by Gilmore only for 
three initial dates, has not been taken into account in the fitting 
procedure as its inclusion degrads the fit quality. This is not 
suprising as the U band is dominated by the Balmer continuum 
which is very sensitive to nonstandard phenomena such as de- 
partures from hydrostatic equilibrium, non-LTE and winds. As 
can be seen from Fig. (21 V4332 Sgr was systematically brighter 
by ~0.5 magnitude in the U band when compared to the stan- 
dard supergiants and giants. 

Note also that when fitting the magnitudes derived from 
IMWT99l for 5-6 June, the B magnitude was not taken into ac- 
count. It was presumably affected by strong emission lines, as 
can be seen from Fig. 5 in ' MWT99I The results obtained for 
this date are particularly uncertain and should be regarded with 
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Fig. 3. The specti-um of V4332 Sgr in 1998, 1999 and 2003. 
Asterisks - JHK in May 1998 from 2MASS. Triangles 
- IcunnJK, in September 1999 from DENIS. Circles - 
BVRJcJHK on May-September 2003 (see the text). Full and 
dotted curves - standard supergiant and giant (respectively) 
spectra fitted to the 1999 and 2003 data (for details see the text). 
Parameters of the fits are given in Tabled The model spectra 
have been reddened with Eb-v - 0.32. 



caution. They are based on rough estimates of the VRI mag- 
nitudes made from the spectral observations. Besides, judging 
from the estimated effective temperature the observations cov- 
ered only the shortwavelength Wien's part of the spectrum. The 
bulk of the energy of the object was presumably emitted in the 
infrared where no measurements were made. Finally, our fits 
for this date result from extrapolation beyond the range of the 
standard spectra (the latest types for which intrinsic colours are 
available are: M6-7 for giants and M5-6 for supergiants). 

Figure |3] shows the results of observations in 1998, 1999 
and 2003. Asterisks indicate the 2MASS measurements ob- 
tained on 18 May 1998. Triangles show the DENIS results de- 
rived on 11 September 1999. Circles represent the spectrum 
observed in May-September 2003 (see above for the sources fo 
the data). 

First we discuss the observations obtained in 2003 as they 
cover the largest spectral range. The whole spectrum, i.e. all the 
circles in Fig. |3] cannot be fitted with a single standard spec- 
trum. This can only be done for shorter wavelengths. In Fig.|3 
we show the best fits of the supergiant (full curve) and giant 
(dotted curve) spectra for the VRch-J measurements. The pa- 
rameters of the fits are given in Table |3] The B magnitude has 
not been taken into account in the fitting procedure because of 
its significant uncertainty. However, as can be seen from Fig.|3 
it fits well the obtained spectra. In the long wavlength range the 
H and, particularly, K magnitudes show a clear excess com- 
pared to the spectra fitted in the shorter wavelengths. With the 
L and M magnitudes measured bv .BVA04. (not shown in Fig.|3j 



one can easily conclude that the source of this excess dominates 
the brightness of the object in the infrared. This infrared excess 
will be discussed in Sect.|6l 

The infrared observations displayed in Fig. |3l show that 
V4332 Sgr evolved systematically between 1998 and 2003. It 
was becoming fainter in J (and /) but brighter in K. This can be 
interpreted as due to a gradual fading of the main object (seen 
in J and /) and the increasing infrared excess in K. Therefore 
when fitting the standard spectra to the DENIS magnitudes ob- 
tained on 11 September 1999 (triangles in Fig. 0, expecting 
that K may be affected by the IR excess, we considered only / 
and J. Indeed, as can be seen from Fig.|3] the K magnitude is 
~0.5 magnitude above the fitted spectrum. Note, however, that 
if all the three (UK) measurements are taken into account the 
fits would be quite poor and Tetj would lower by ~100K for 
class 111 and ~200 K for class 1 compared to the values given in 
TableE 

No fit to the 2MASS data on 18 May 1998 (asterisks) is 
shown in Fig. |3] The spectral range of the data is quite nar- 
row and hence any Teff between ~3000 K and -4500 K could 
be considered to satisfactorily fit the data. Fortunately for this 
wide Teff range, log 9 of the fits varies by less than 0.20. The 
mean values of log 6 from these fits are given in Table|3] 

4.4. Evolution of the object 

The last two columns in Table|3lshow the stellar radius and lu- 
minosity (given in solar units) of V4332 Sgr calculated from 
the angular radius (6) and the eff'ective temperature (Teff) as- 
suming a distance of 1.8 kpc. 

Figure 13 displays the evolution of principal parameters of 
V4332 Sgr taken from Tabled The time is in days since 24 
February 1994, i.e. since the discovery of V4332 Sgr in erup- 
tion. The effective temperature is in units of 10^ K, while the 
stellar radius and luminosity are in solar units. Open symbols 
show the results from fitting the standard supergiant (1) spectra 
to the observations. The same but using the giant (111) spectra is 
presented with the full points. Asterisks in panel (a) denote the 
effective temperatures derived by MWT99 from fitting stellar 
atmosphere models to their spectra. 

Let us first discuss the spectral types obtained from our 
analysis. As can be seen from Table|3]no matter which luminos- 
ity class of the standard spectra (I or III) is used, the resultant 
spectral types are always very similar. Moreover, our spectral 
classes are also very close to those obtained by MWT99 from 
the classification of their spectra. This shows that the general 
approach adopted in our analysis of the photometric data is 
consistent and reliable. 

The consistency of the effective temperatures derived in the 
different ways is not as good as that of the spectral types. As 
can be seen from Tableland Fig.l^aj initially, when the object 
was of the K type, both kinds of standard spectra give practi- 
cally the same values of Teff. Later, when the object was be- 
coming cooler, an increasing discrepancy appears in the sense 
that Teff obtained using the supergiant spectra is systematically 
lower than that from using the giant spectra. The reason lies in 
the cahbrations of Jeff versus spectral type used in our study 
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Fig. 4. Evolution of V4332 Sgr with time. Time (displayed in log scale) is in days counted since the discovery of the object in 
eruption (24 Feb. 1994). Left (a) - effective temperature in lO-' K. Centre (b) - logarithm of radius in Rq. Right (c) - logarithm 
of luminosity in Lq. Open symbols - results from fitting the supergiant (I) spectra. Full points - results from fitting the giant (III) 
spectra. Asterisks in (a) - effective temperature from MWT99 Asterisks in (b) and (c) - estimates on 5-6 June 1994 based on 
the results of lMWT^^ as discribed in the text. Full curves in (b) and (c) - results of a model discribed in Sect.|5] 



Table 3. Evolution of V4332 Sgr in 1994-2003 

Date Sp.type -log R/Rq L/Lq 



(Schmidt-Kaler [T982t . For the same spectral type, the value of 
Tetf for supergiants is systematically lower than that for the gi- 
ants but for the G, K types the difference is small. For the M 
types the discrepancy increases and becomes as large as 530 K 
at M5. This has obvious consequences for the derived effec- 
tive radii and luminosities: the use of the supergiant standard 
spectra usually results in larger radii and luminosities. 



Comparing our results with those obtained bv lMWT99l one 
sees that on March 4 and 9 their T^ff is systematically higher. 
This agrees with what was noted in IMWT991 namely that the 
value of Tetf they obtained from the stellar model analysis was 
too high for the observed spectral type. MWT99 ascribed this 
discrepancy to the narrow observed spectral range. On the later 
dates the spectral types and T^ff derived in iMWT99> are consis- 
tent and, as can be seen from Fig. their T^g follows quite 
closely our values obtained using the supergiant spectra. 

Note that, as discussed in Sect. 14.31 our results obtained 
from the data on 5-6 June 1994 (symbols at log t ^ 2.0 in 
Fig. |4} are subject to significant uncertainty. They are based 
on uncertain observational estimates done only in the Wien's 
part of the spectrum, as well as, on the extrapolated standard 
spectra. Due to the increasing discrepancy in the effective tem- 
perature between giants and supergiants for the late M types, 
discussed above, the difference between our estimates of Teff is 
as large as 900 K. This resulted in large differences in the stel- 
lar radius and luminosity seen in Table|3land Fig.|3 MW T99I 
from their model atmosphere analysis, estimated that by 5- 
6 June the object had faded by a factor of 100 in luminosity 
compared to the beginning of March. Adopting our determi- 
nations of ~5000Lo for '^he beginning of March this results 
in an estimate of ~5OL0 and, assuming Teff - 2300 K (from 
IMWT99> . an effective radius of -44 Rq on 5-6 June. These 
estimates are shown with asterisks in Fig.0Jfe)fc). 

As can be seen from Table|3]and Fig.|3; during ~10 days 
after its discovery, V4332 Sgr sustained a fairly constant lumi- 
nosity, the star being ~5000 times brighter than before the erup- 
tion. The radius was then ~ 140 times larger than that of the pro- 
genitor, possibly slowly expanding with time. Near 5 March, 
the stellar luminosity declined rapidly and during ~6 days the 
object faded by a factor ~1. 

Concerning the effective radius we can only state, given the 
uncertainties in our results, that the radius was not changing 
significantly till 10 March. It seems, however, that a significant 



27.7Feb.'94 K0.2 1 4400. 8.85 114. 4360. 

K1.7III 4470. 8.86 110. 4380. 

3.7Mar.'94 K3.3 1 4050. 8.75 142. 4880. 

K3.6III 4070. 8.75 141. 4950. 

4.5 Mar.' 94 K3-4 III-I 4400. M. 

4.7Mar'94 K4.0 1 3950. 8.73 150. 4910. 

K4.2III 3970. 8.73 149. 4980. 

7.7Mar'94 Ml. II 3640. 8.75 141. 2790. 

Ml. 1 III 3710. 8.80 128. 2790. 

8.7Mar'94 M2.4 1 3340. 8.71 154. 2690. 

M2.7III 3560. 8.78 133. 2550. 

9.5Mar'94 M3 III-I 3800. M. 

10.7Mar'94 M3.6 1 3070. 8.67 169. 2290. 

M4.1III 3420. 8.82 119. 1760. 

11.4Mar'94 M5-6 III-I 3100. M. 

M3.8I 3040. 8.89 102. 790. 

M4.1III 3420. 9.08 66. 540. 

20.4Mar'94 M6-7 III-II 2600. M. 

5-6Jun.'94 M8-9 III 2300. M. 

M9.0I 2120. 8.97 86. 135. 

M8.0III 3050. 9.87 10.7 8.9 

18 May '98 9.96 8.7 

10.06 6.9 

11 Sep.'99 M3.9I 3000. 9.88 10.6 8.2 

M4.4III 3390. 10.06 6.9 5.6 

May-Sept.'03 M2.7 1 3280. 10.13 5.9 3.6 

M2.9III 3540. 10.22 4.8 3.2 
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shrinkage of the radius occurred near 1 1 March. The initial lu- 
minosity decline was thus primarily due to the decline in the 
effective temperature. 

The subsequent fading involved the decline of all the stellar 
parameters, i.e. L, Teff and R. At a certain time T^ff reached a 
minimum but the available data do not permit an estimate of 
when it occurred and how deep this minimum was. From Fig.0] 
it can be inferred that this occured near June 1994 or later. The 
lowest observed T^f[ was reached on 5-6 June 1994, i.e. some 
90 days after the luminosity started declining. During this time 
interval, the object dropped in luminosity by factor at least ~40. 
According to MWT99 this luminosity drop was of factor 100 
and the object became as cool as 2300 K. 

The decline in luminosity and radius has continued, al- 
though at a decreasing rate. This has been followed by a slow 
rise in the effective temperature and at present the object is sig- 
nificantly hotter than in June 1994. Comparing the results ob- 
tained in 2003 with the estimates made for the progenitor (see 
Sect. 14. 2> . we can conclude that V4332 Sgr is now larger and 
more luminous, but cooler than before the eruption, so it has 
not yet returned to its pre-outburst state. 

5. Gravitational contraction in the decline 

The fading of V4332 Sgr can be consistently interpreted if one 
assumes that the energy source, that caused V4332 Sgr to erupt 
to ~140 Rq, suddenly disappeared. This resulted in a gravita- 
tional contraction of the inflated envelope and a release of its 
gravitational energy. To show this, let us assume that the in- 
flated envelope has a density distribution of p ~ r^^^^ (follow- 
ing Soker & Tylenda 2003 1. Then, assuming that the envelope 
mass. Me, is small compared to the star mass, M*, it is straight- 
foward to show that the gravitational energy of the envelope is 



iR^Ry^' 



(1) 



where G is the gravitational constant, Ri, is the radius of the 
unperturbed stellar layers (assumed to be the inner radius of 
the envelope) and R is the outer radius of the envelope. If the 
radiated luminosity, L, is due to gravitational contraction of the 
envelope, one gets 



L = - 



dEg _ GM^M^ 1 dR 
"dT ~ 2rI!^ Wi^~dt' 

On the other hand 



L = AnR^o-T^^. 



Combining Eq. ^ with Eq. gives 



dR 8KR'f(rT^^^. 



RV2 



GM^M^ 



dt. 



(2) 



(3) 



(4) 



Assuming that Teff is constant (which is obviously a very crude 
approximation), Eq. can be integrated giving 



1 _ lOnR'fcrT^g 1 



(5) 



where Rq is the initial envelope radius assumed to start con- 
tracting at f = 0. 

The solid curves in Fig.l^fejfc) shows the evolution of a 
gravitationally contracting envelope calculated from Eq. (|5}- 
The solar values for the stellar paramaters, i.e. M*. = I.QMq 
and Ri, - I.QRq, have been assumed, as well as, Tetf = 3300 K 
as a typical effective temperature during the decline. The mass 
of the envelope have been chosen to fit the observed evolution 
and in the case shown in ¥\g.^b)(c), is 1.5 x \0'^Mq. 
We note that values of Mg differing by less than a factor 3, i.e. 
being within 0.5 - 5.0 x IO^^Mq, can be considered as reason- 
ably reproducing the observations. Given the crudness of our 
approach and uncertainties in the observational determinations 
of R, we can conclude that pure gravitational contraction of the 
envelope inflated during the eruption can satisfatorily explain 
the observed decline of V4332 Sgr The mass of the envelope 
is small, of order IO^^Mq. Note that the total mass involved in 
the eruption was probably significantly larger than Mg. Some 
mass was probably lost from the object and/or is circulating 
the object in a disc-like structure (see Sect.|6}. 

6. Infrared excess in 2003 

One of the most important features in the spectrum of 
V4332 Sgr in 2003 is the infrared excess clearly seen in the 
H and K bands in Fig. |3] As discussed in Sect. 14.31 no such 
excess is necessary to explain the data in 1998. The 2MASS 
magnitudes can be well accounted for by a cool star some- 
what brighter then the stellar component seen in 2003. In 
September 1999 the K band was probably affected by an ex- 
cess, although its magnitude was certainly significantly smaller 
than that in 2003. 

In order to discuss the origin of this IR excess, we 
have attempted to fit the observations in 2003 with multi- 
component spectra. Symbols in Fig. |5] show the spectrum in 
May-September 2003 (see Sect. l4.3l for the sources of the data). 
The full curve shows a fit of the data with a two component 
spectrum: a standard supergiant spectrum of T^s - 3280 K 
(M2.7 I - the same as in Fig.|3} and a blackbody distribution 
with Tbb = 750 K. As can be seen from Fig.|5lthe fit is very 
good. The angular radius of the blackbody component in Fig.|3] 
is log 9 - -8.28, which at a distance of 1.8 kpc gives an effec- 
tive radius of ~AQQRq and a luminosity of ~50 Lq. 

Fits of the kind shown in Fig.|5lcan be interpreted as evi- 
dence for a dusty envelope that reprocesses radiation absorbed 
from the central star. This interpretation is however very un- 
likely. The infrared component is ~15 times brighter than our 
estimate of the stellar luminosity in Table |3l This would im- 
ply that the central star is significantly more luminous and that 
what we observe in the optical is the stellar spectrum severely 
attenuated by dust. Then the central star, reddened by dust, 
would be much hotter than estimated from the photometry. Yet, 
as noted in Sect.|3] our spectrum suggests an early M type star, 
thus consistent with the photometric results. Besides, the large 
present luminosity and the presence of hot nearby dust would 
inevitably indicate that V4332 Sgr has recently experienced 
another outburst. With the effective radius of the blackbody 
component estimated above (~4OO^0) and adopting a rather 
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Fig. 5. The magnitudes observed in May-September 2003 (see 
Sect. 14. 31 for the sources of the data) fitted with a two compo- 
nent spectrum: M2.7 supergiant of Tetf - 3280 K and black- 
body of r = 750 K. Symbols - observed magnitudes. Dotted 
curves - individual component spectra. Full curve - sum of the 
two components. 



lower Umit of 10 km/s for the expansion velocity, one finds that 
dust would have been lost a year before May-September 2003. 
However, as discussed in Sect. 14.31 the infrared excess started 
developing in 1999. Besides, there is no other observational in- 
dication that the object recently suffered from a secondary out- 
burst. On the contrary, as can be seen from Fig.|3 the object has 
been following the long term decline after the 1994 eruption. 
The lack of any significant evolution in the optical magnitudes 
between May and September 2003 is also an evidence for the 
long time scale evolution of the object in 2003. 

BVA04 have recently considered that V4332 Sgr might be 
a young object surrounded by a protostellar circumstellar disc. 
This hypothesis has been invoked to explain the origin of the 
water-ice band observed in the infrared spectrum. It also allows 
[BVA04 to suggest, following the ideas of Soker & Tylenda 
( l2003l and Retter & Marom (l2003l . that the 1994 outburst of 
V4332 Sgr might have been due to infall of an inner planet 
due to tidal interaction with the disc. According to BVA04 the 
presently observed emission spectrum and the infrared excess, 
which appeared in 1999, might have been due to interaction of 
the 1994 ejecta with the circumstellar disc. However, an inter- 
action of this kind should have happen during the 1994 eruption 
or soon after it. Assuming an expansion velocity of 100 km/s 
( MWT99l the distance of 400 Rq is reached in a month. Thus 
it is not clear why the IR excess might have appeared ~5 years 
later and why it might have been brightening with time, as ob- 
served. 

The idea of a protostellar disc can however be used to ex- 
plain the increasing IR excess, although in a different way to 
that in IBVA04i In the case of an inner planet and an external 



accretion disc, tidal forces transfer angular momentum from 
the orbital motion of the planet to the inner edge of the disc. As 
a result the inner edge of the disc is kept at a certain distance 
from the planet orbit, determined by a balance between the rate 
of the angular momentum gain from the planet and the rate 
in which the angular momentum is transported outward in the 
disc. When the planet gets accreted by the central star the erup- 
tion can disturb the disc, possibly destroying its inner regions. 
After the event the inner disc does not gain angular momen- 
tum any more. As a result the inner edge will be approaching 
the star and the inner disc will become hotter and hotter. Thus 
we can speculate that just after the 1994 eruption the inner disc 
in the V4332 Sgr system was far away, cold and thus radiat- 
ing only in the far IR, while at present the inner disc regions 
are closer to the star, thus they are hotter and dominating the 
observed brightness in the KLM bands. 

There are, however, observational facts which do not favour 
the hypothesis that V4332 Sgr is a young object. First, the 
position of V4332 Sgr in galactic coordinates is / - 13°63, 
b = -9?40. According to catalogues of Sharpless il959> . Lynds 
(jT965^, Avedisova (2002i and Russeil (2003) there is no HII 
or star forming region closer than 5° from the position of 
V4332 Sgr. Also from the CO map of Dame et al. ( dooTI Fig 2) 
one can see that the object lies well outside any significant CO 
emission. Second, .MWT99. from their high resolution spec- 
tra near Hq,, have derived a radial velocity of V4332 Sgr to 
be -180 km/s. A similar value, i.e, -160 km/s, can be ob- 
tained from the observed positions of all the emission lines 
listed in their Table 3. This result has been derived from out- 
burst spectra so interpretations other than the radial velocity 
of the object (although very unlikely) can be considered (see 
IMWT99l l. However, the observed wavelengths of the Cal line 
in our Table|2still gives a similar radial velocity, i.e - 140 km/s. 
The Galactic rotation curve (see e.g. Brand & Blitz . 1993j pre- 
dicts, for the position of V4332 Sgr and a distance of 1.8 kpc, 
Vlsr = -1-13 km/s which is equivalent to a heliocentric radial 
velocity of ~ -1-3 km/s. A CO radial velocity map of Dame at al. 
|2001 Fig. 3) shows, at the galactic longitude of V4332 Sgr, 
Vlsr between and +140 km/s. Thus V4332 Sgr does not fol- 
low the Galactic rotation which is not expected to be the case 
for a young object. 

Formation of a circumstellar disc is, however, also possible 
in the merger scenario proposed by Soker & Tylenda (2003 1 
in which a star accretes a less massive object. In this case, 
the inflated stellar envelope gains not only energy from the 
merger, but also angular momentum from the orbital motion. 
Note that if the merger event is provoked by an interaction with 
a fast moving object (as discussed above V4332 Sgr is prob- 
ably moving with a velocity ^200 km/s relative to the objects 
rotating in the Galactic disc) the angular momentum accreted 
is likely to be much larger than that of the Keplerian motion 
near the star surface. During the outburst, i.e. when the enve- 
lope remains inflated to large radii, processes due to turbulent 
or convective motions, tidal effects, as well as magnetic fields 
are likely to transport angular momentum outward. When the 
merger process is more or less complete, the inflated envelope 
would tend to collapse but the angular momentum stored in the 
equatorial regions may prevent these regions from contracting 
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significantly. They would form a ring-like or disc-like structure 
orbiting the central object. During the fading phase when the 
contracting central object is likely to rotate significantly faster 
than the matter orbiting at larger distances, magnetic fields and 
tidal forces can continue to transport angular momentum out- 
ward. With time, due to viscous processes, the ring-like struc- 
ture would evolve toward an accretion disc and dissipate its 
energy. 

Following the above discussion we have fitted the observa- 
tions in 2003 with a composite spectrum of a central star and a 
standard blackbody accretion disc. The central star spectrum is 
the same as in Fig.|5l The distribution of the effective tempera- 
ture of the disc is (see e.g. Eq. 5.43 in Frank et al. l2002> 



T(R)^To{3 



where 



1/4 



7^0 = 



( GM,M' 



1/4 



(6) 



(7) 



/?() is the inner radius of the disc while M is the accretion 
rate. Note that T{R) attains a maximum value 0.642 To at 
R - 1.36 /?o- Integrating the blackbody spectrum with T(R) 
over the disc surface one gets the disc spectrum. 

We do not show the final fit to the observations as it is 
very much the same as that in Fig. |5] It has been obtained 
for To = 1400 K and log(6lo) = -8.70, where 0o is the "ob- 
served" angular inner radius of the disc. Assuming a distance 
of 1.8 kpc and M, = 1 Mq this translates to Rq ^ 160 T?© 
and M ^ 8.8 X lO^'^Mo/yr. The total luminosity of the disc is 
then ~87 Lq. Note that is very close to the maximum radius 
of the star during the 1994 eruption (see Table |3}- This sug- 
gests that if the disc is of protostellar origin its inner regions 
have probably been distroyed during the eruption, as discussed 
above. 

Taking Eq (5.69) in Frank et al. 1120021 and assuming the 
viscosity parameter, a - 0.01 - 0.1, one gets, for the above 
values of R^ and M, the viscous time scale of 4-20 years. Thus 
one can argue that in 1998-1999 the inner edge of the disc or 
ring-like structure was somewhat farther away and cooler than 
in 2003, the star was brighter so the relative contribution of the 
structure in the HK bands was negligible. 

7. Discussion and conclusions 

There is no doubt that the eruption of V4332 Sgr observed in 
1994 was a highly unusual, enigmatic and extremely interest- 
ing astrophysical event. This statement results not only from 
the observed behaviour of the object in 1994 but also from 
its present state marked by the unusual emission line spec- 
trum and the mysterious infrared excess. It is most unfortunate 
that so little information is available regarding the evolution 
of V4332 Sgr during and after its 1994 eruption. The data is 
limited to the optical with just the six multicolour photomet- 



was likely emitting much of its energy in this wavelength re- 
gion. 

The results of our photometric analysis in Sect. |4] supple- 
mented with the results of .MWT99. clearly show that the de- 
cline of V4332 Sgr following its 1994 eruption was marked 
by a decline in all principal parameters: luminosity, radius and 
effective temperature. In 2003 the object was ~1500 times 
less luminous than during the eruption but still remaining in 
the spectral range of the M type. As discussed in Soker & 
Tylenda (2 003>) . this type of evolution rules out the classical 
nova mechanism which could have worked for a hypotheti- 
cal white dwarf companion to the progenitor. Note also that in 
June 1994 and May-September 2003, V4332 Sgr was signifi- 
cantly (1.0-1.5 mag.) fainter in the B band than the progenitor 
(between these two dates there were no measurements in E). 
This implies that it was the progenitor which erupted in 1994 
and not its hypothetical faint companion. 

As discussed in Sect. |5] the decline of V4332 Sgr can be 
understood in terms of a pure gravitational contraction of the 
inflated stellar envelope. An event of this kind is expected to 
follow a sudden switch-off of the energy sources which have 
so far been working at the base of the envelope. This interpre- 
tation is well in line with the stellar merger scenario of Soker & 
Tylenda (2003 1. When the merger phase is complete, the source 
of energy quickly drops and the inflated envelope has to con- 
tract under the gravity of the star The considerations in Sect.|5] 
show that in order to explain the observed decline of V4332 Sgr 
over almost a decade, the mass of the inflated stellar envelope 
must be small, of order 1O"^M0. 

The origin of the infrared excess dominating the observed 
KLM magnitudes in 2003 is unclear. When fitted with black- 
body it implies a temperature of ~750 K and a luminosity ~15 
times higher than that of the central star As discussed in Sect.|6l 
it is very unlikely that the excess is due to a dusty envelope 
reprocessing the radiation from the central star We have con- 
sidered that the IR radiation can originate from a circumstellar 
disc dissipating its gravitational energy via viscous processes. 
The disc can be either of protostellar origin, following sug- 
gestions of IB VA04I or being a remnant storing angular mo- 
mentum after the 1994 eruption. Both possibilities are in line 



ric observations of Gilmore J1994> covering 1 1 days and seven 
spectra of IMWT99 ' obtained in a three month period. The lack 
of infrared data is particularly disappointing since the object 



with the merger scenario of Soker & Tylenda J2003> . The first 
one assumes that V4332 Sgr is a young object. This hypoth- 
esis is however inconsistent with the position of the object in 
Galaxy and, particularly, with its radial velocity significantly 
different from the Galactic rotation curve. Future observations 
of V4332 Sgr in the infrared, particularly in longer wavelengths 
might be crucial for understanding the nature of the object. 

A lower limit to the mass of this object accreted in the 
merger scenario of Soker & Tylenda (2Q03) can be obtained 
from the total energy emitted in the eruption. The luminosity of 
V4332 Sgr integrated since the discovery of the 1994 eruption 
till 2003 is ~ 4.5 x 10'*^ ergs. Equating this to GM*Macc//?* 
one obtains (assuming - I.OMq and = l.O/?©) 
A^acc - IO^^Mq. Note that this estimate does not take into 
account the energy radiated away before the discovery (it is 
very likely that the eruption of V4332 Sgr started well before 
28 February 1994), nor the energy stored in the matter now 
circulating the central star, nor the kinetic energy in mass loss. 
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Given the above estimate as well as that of Me in Sect.|5] we can 
conclude that the mass of the accreted object was ^ 10""* Mq. 

The unusual emission line spectrum observed in April 2003 
and described in Sect. 13 deserves futher studies as it can pro- 
vide important insight into the current state of V4332 Sgr From 
our preliminary analysis in Sect.|3we can conclude that it must 
originate in an optically thin, neutral, molecular, cold medium. 
Low optical thickness means low column density which, with 
large equivalent widths of many lines in the spectrum, requires 
large volume - considerably larger than that of the central 
star. This obviously rules out the stellar atmosphere and re- 
gions in its near vicinity. The low rotational temperature in- 
ferred in Sect. |3l (from different emission bands) implies that 
the medium, in the bulk, is significantly cooler than lOOOK. 
This is the sort of temperature we have found for the disc-like 
structure in Sect. |6] However, this structure is expected to be 
mostly optically thick. It seems that possible sites for the emi- 
sion line spectrum to originate might be in regions between 
the stellar surface and the disc where the contracting stellar en- 
velope might have left some matter circulating now near the 
equatorial plane. IBA04I have found that the observed width of 
the KI emission lines if interpreted with Doppler broadening 
give a velocity of 260 km/s. The measured FWHM of the Cal 
line in our spectrum (see Table |2jl when corrected for the in- 
strumental FWHM (5.6 A) gives a similar value, i.e. 280 km/s. 
This value would correspond to a Keplerian velocity at ~10 Rq 
for a 1 Mq star The required excitation could result from the 
stellar radiation or viscous processes, or both. 

The galactic position of V4332 Sgr, I = 13?63 b = -9?40, 
might suggest that the object is related to the Galactic bulge. 
Thus V4332 Sgr would be at a distance of ~ 8.5 kpc, which, 
on the one hand, seems to be too large for the observed reden- 
ning, as discussed in Sect. 14.21 On the other hand, however, it 
would be easier in this case to understand the observed radial 
velocity of V4332 Sgr, which is, as discussed in Sect. |6j in- 
consistent with the rotation of the Galactic disc. At a distance 
of 8.5 kpc all the radii and luminosities derived in this paper 
would increase by a factor of 4.7 and 22, respectively. In par- 
ticular, the progenitor would be a G type subgiant of a mass 
of ~ 2 Mq, as can be inferred from stellar evolutionary tracks 
(e.g. Iben ll965l . evolving from the main sequence to the red 
giant branch. The star would have a ~ 0.2 Mq helium core 
and an envelope with a mean density of 0.02 g cm"^. Thus one 
could refer to the scenario of Retter & Marom ( 2003 ) while in- 
terpreting the eruption of V4332 Sgr. Indeed Mace - 1O"'*M0 
is then required to explain the energy emitted in the eruption 
so a massive planet or a brown dwarf would be involved in the 
merger event. However, the accreted object, being significantly 
denser than the subgiant envelope, would be expected to pene- 
trate deeply in the envelope, as discussed e.g. in Livio & Soker 
J1984> . Thus most of the subgiant envelope would be disturbed 
and the whole event, the contraction phase in particular, would 
precede on a much longer time scale than observed. Indeed, 
modelling analogous to that in Sect. |5l but done to match the 
observed increase of the radius by a factor of 4.7, assuming 
a central star mass of 2 Mq, requires Mg ^ 3 x 10^'^Mq. In 
the discussed case we would rather expect Me - IM0 to be 
involved. Note that Livio & Soker (il984< show that a merger 



event in the case of a giant and a massive planet may last sev- 
eral thousands years. 

There is no doubt that futher detailed observations of 
V4332 Sgr are necessary. Photometric and spectroscopic mea- 
surements in the near and far infrared would be of particular 
value as the object is emitting most of its energy in this wave- 
length range. High quality spectroscopy in the optical would 
also be very important for understanding the origin of the un- 
usual emission-line spectrum as well as for analysing the un- 
derlying stellar component. In particular, high resolution spec- 
troscopy in the atomic lines would allow precise determination 
of the radial velocity of the object, as well as investigation of 
the kinematic structure of the emitting region. The object is 
probably evolving on a time scale of years so long term ob- 
servational monitoring should be undertaken. We hope that the 
analysis of the existing observational data in this paper, fol- 
lowed by our interpretation and conclusions (which in several 
points may be regarded as speculative) will stimulate deeper 
astrophysical interest in V4332 Sgr, V838 Mon and, possibly, 
other objects of similar nature. 
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